Antagonistic control of rice immunity against distinct pathogens by the two transcription modules via salicylic acid and jasmonic acid pathways.

Although the antagonistic effects of host resistance against biotrophic and necrotrophic pathogens have been documented in various plants, the underlying mechanisms are unknown. Here, we investigated the antagonistic resistance mediated by the transcription factor ETHYLENE-INSENSITIVE3-LIKE 3 (OsEIL3) in rice. The Oseil3 mutant confers enhanced resistance to the necrotroph Rhizoctonia solani but greater susceptibility to the hemibiotroph Magnaporthe oryzae and biotroph Xanthomonas oryzae pv. oryzae. OsEIL3 directly activates OsERF040 transcription while repressing OsWRKY28 transcription. The infection of R. solani and M. oryzae or Xoo influences the extent of binding of OsEIL3 to OsWRKY28 and OsERF040 promoters, resulting in the repression or activation of both salicylic acid (SA)- and jasmonic acid (JA)-dependent pathways and enhanced susceptibility or resistance, respectively. These results demonstrate that the distinct effects of plant immunity to different pathogen types are determined by two transcription factor modules that control transcriptional reprogramming and the SA and JA pathways.

Rice-Magnaporthe oryzae interactions in resistant and susceptible rice cultivars under panicle blast infection based on defense-related enzyme activities and metabolomics.

Rice blast, caused by rice blast fungus (Magnaporthe oryzae), is a global threat to food security, with up to 50% yield losses. Panicle blast is a severe form of rice blast, and disease responses vary between cultivars with different genotypes. Reactive oxygen species (ROS)-mediated signaling reactions and the phenylpropanoid pathway are important defense mechanisms involved in recognizing and resisting against fungal infection. To understand rice-M. oryzae interactions in resistant and susceptible cultivars, we determined dynamic changes in the activities of five defense-related enzymes in resistant cultivar jingsui 18 and susceptible cultivar jinyuan 899 infected with M. oryzae from 4 to 25 days after infection. We then performed untargeted metabolomics analyses to profile the metabolomes of the cultivars under infected and non-infected conditions. Dynamic changes in the activities of five defense-related enzymes were closely related to panicle blast resistance in rice. Metabolome data analysis identified 634 differentially accumulated metabolites (DAMs) between resistant and susceptible cultivars following infection, potentially explaining differences in disease response between varieties. The most enriched DAMs were associated with lipids and lipid-like molecules, phenylpropanoids and polyketides, organoheterocyclic compounds, organic acids and derivatives, and lignans, neolignans, and related compounds. Multiple metabolic pathways are involved in resistance to panicle blast in rice, including biosynthesis of other secondary metabolites, amino acid metabolism, lipid metabolism, phenylpropanoid biosynthesis, arachidonic acid metabolism, arginine biosynthesis, tyrosine metabolism, tryptophan metabolism, tyrosine and tryptophan biosynthesis, lysine biosynthesis, and oxidative phosphorylation.

Comparison of Transcriptome between Tolerant and Susceptible Rice Cultivar Reveals Positive and Negative Regulators of Response to Rhizoctonia solani in Rice.

Rice (Oryza sativa L.) is one of the world's most crucial food crops, as it currently supports more than half of the world's population. However, the presence of sheath blight (SB) caused by Rhizoctonia solani has become a significant issue for rice agriculture. This disease is responsible for causing severe yield losses each year and is a threat to global food security. The breeding of SB-resistant rice varieties requires a thorough understanding of the molecular mechanisms involved and the exploration of immune genes in rice. To this end, we conducted a screening of rice cultivars for resistance to SB and compared the transcriptome based on RNA-seq between the most tolerant and susceptible cultivars. Our study revealed significant transcriptomic differences between the tolerant cultivar ZhengDao 22 (ZD) and the most susceptible cultivar XinZhi No.1 (XZ) in response to R. solani invasion. Specifically, the tolerant cultivar showed 7066 differentially expressed genes (DEGs), while the susceptible cultivar showed only 60 DEGs. In further analysis, we observed clear differences in gene category between up- and down-regulated expression of genes (uDEGs and dDEGs) based on Gene Ontology (GO) classes in response to infection in the tolerant cultivar ZD, and then identified uDEGs related to cell surface pattern recognition receptors, the Ca2+ ion signaling pathway, and the Mitogen-Activated Protein Kinase (MAPK) cascade that play a positive role against R. solani. In addition, DEGs of the jasmonic acid and ethylene signaling pathways were mainly positively regulated, whereas DEGs of the auxin signaling pathway were mainly negatively regulated. Transcription factors were involved in the immune response as either positive or negative regulators of the response to this pathogen. Furthermore, our results showed that chloroplasts play a crucial role and that reduced photosynthetic capacity is a critical feature of this response. The results of this research have important implications for better characterization of the molecular mechanism of SB resistance and for the development of resistant cultivars through molecular breeding methods.

Genome-Wide Association Study Identifies a Plant-Height-Associated Gene OsPG3 in a Population of Commercial Rice Varieties.

Plant height is one of the most crucial components of plant structure. However, due to its complexity, the genetic architecture of rice plant height has not been fully elucidated. In this study, we performed a genome-wide association study (GWAS) to determine rice plant height using 178 commercial rice varieties and identified 37 loci associated with rice plant height (LAPH). Among these loci, in LAPH2, we identified a polygalacturonase gene, OsPG3, which was genetically and functionally associated with rice plant height. The rice plant exhibits a super dwarf phenotype when the knockout of the OsPG3 gene occurs via CRISPR-Cas9 gene-editing technology. RNA-Seq analysis indicated that OsPG3 modulates the expression of genes involved in phytohormone metabolism and cell-wall-biosynthesis pathways. Our findings suggest that OsPG3 plays a vital role in controlling rice plant height by regulating cell wall biosynthesis. Given that rice architecture is one of the most critical phenotypes in rice breeding, OsPG3 has potential in rice's molecular design breeding toward an ideal plant height.

Identification of two novel rice S genes through combination of association and transcription analyses with gene-editing technology.

Traditional rice blast resistance breeding largely depends on utilizing typical resistance (R) genes. However, the lack of durable R genes has prompted rice breeders to find new resistance resources. Susceptibility (S) genes are potential new targets for resistance genetic engineering using genome-editing technologies, but identifying them is still challenging. Here, through the integration of genome-wide association study (GWAS) and transcriptional analysis, we identified two genes, RNG1 and RNG3, whose polymorphisms in 3'-untranslated regions (3'-UTR) affected their expression variations. These polymorphisms could serve as molecular markers to identify rice blast-resistant accessions. Editing the 3'-UTRs using CRISPR/Cas9 technology affected the expression levels of two genes, which were positively associated with rice blast susceptibility. Knocking out either RNG1 or RNG3 in rice enhanced the rice blast and bacterial blight resistance, without impacting critical agronomic traits. RNG1 and RNG3 have two major genotypes in diverse rice germplasms. The frequency of the resistance genotype of these two genes significantly increased from landrace rice to modern cultivars. The obvious selective sweep flanking RNG3 suggested it has been artificially selected in modern rice breeding. These results provide new targets for S gene identification and open avenues for developing novel rice blast-resistant materials.

PacBio full-length transcriptome analysis provides new insights into transcription of chloroplast genomes.

Chloroplast and mitochondrial DNA (cpDNA and mtDNA) are apart from nuclear DNA (nuDNA) in a eukaryotic cell. The transcription system of chloroplasts differs from those of mitochondria and eukaryotes. In contrast to nuDNA and animal mtDNA, the transcription of cpDNA is still not well understood, primarily due to the unresolved identification of transcription initiation sites (TISs) and transcription termination sites (TTSs) on the genome scale. In the present study, we characterized the transcription of chloroplast (cp) genes with greater accuracy and comprehensive information using PacBio full-length transcriptome data from Arabidopsis thaliana. The major findings included the discovery of four types of artifacts, the validation and correction of cp gene annotations, the exact identification of TISs that start with G, and the discovery of polyA-like sites as TTSs. Notably, we proposed a new model to explain cp transcription initiation and termination at the whole-genome level. Four types of artifacts, degraded RNAs and splicing intermediates deserve the attention from researchers working with PacBio full-length transcriptome data, as these contaminant sequences can lead to incorrect downstream analysis. Cp transcription initiates at multiple promoters and terminates at polyA-like sites. Our study provides new insights into cp transcription and new clues to study the evolution of promoters, TISs, TTSs and polyA tails of eukaryotic genes.

Identification of a Rice Leaf Width Gene Narrow Leaf 22 (NAL22) through Genome-Wide Association Study and Gene Editing Technology.

Rice leaf width (RLW) is a crucial determinant of photosynthetic area. Despite the discovery of several genes controlling RLW, the underlying genetic architecture remains unclear. In order to better understand RLW, this study conducted a genome-wide association study (GWAS) on 351 accessions from the rice diversity population II (RDP-II). The results revealed 12 loci associated with leaf width (LALW). In LALW4, we identified one gene, Narrow Leaf 22 (NAL22), whose polymorphisms and expression levels were associated with RLW variation. Knocking out this gene in Zhonghua11, using CRISPR/Cas9 gene editing technology, resulted in a short and narrow leaf phenotype. However, seed width remained unchanged. Additionally, we discovered that the vein width and expression levels of genes associated with cell division were suppressed in nal22 mutants. Gibberellin (GA) was also found to negatively regulate NAL22 expression and impact RLW. In summary, we dissected the genetic architecture of RLW and identified a gene, NAL22, which provides new loci for further RLW studies and a target gene for leaf shape design in modern rice breeding.

Full-length chloroplast genome of Dongxiang wild rice reveals small single-copy region switching.

Background: Plant chloroplast DNA (cpDNA) typically has a circular structure, including a large single-copy region (LSC), a small single-copy region (SSC) and two inverted repeats (IR1 and IR2). The organization of these four elementary regions LSC-IR1-SSC-IR2 is highly conserved across all plant cpDNAs. Very few structural variations (SVs) occurring at the elementary-region level have been reported. Results: In the present study, we assembled the full-length cpDNA of Dongxiang wild rice line 159 (DXWR159). Using the long PacBio subreads, we discovered a large inversion of SSC and a large duplication of IR in DXWR159 cpDNAs. Significantly, we reported for the first time forward and reverse SSCs of cpDNAs in similar proportions and named the frequent inversion of a whole SSC as SSC switching. Conclusions: Our study helps researchers to correctly assemble the chloroplast genomes. Our recombination model explained the formation of large SVs in cpDNAs and provided insights into a novel scientific question that if there are common mechanisms in the formation or translocation of all kinds of transposon-like elements (TLEs). We propose that: (1) large inversion is the most accepted mutation type of SVs in cpDNAs; (2) SSC switching ubiquitous occurs in plant cpDNAs; and (3) further investigation of molecular mechanism underlying SSC switching may reveal new driving forces for large SVs.

Genome-edited ATP BINDING CASSETTE B1 transporter SD8 knockouts show optimized rice architecture without yield penalty.

This study reports the identification of the rice open reading frame Semi-Dwarf in chr8 (SD8) that encodes a putative ortholog of Arabidopsis thaliana ABCB1. Genome editing of SD8 leads to optimized rice architecture by reducing plant height and flag-leaf angle without yield penalty. Rice SD8 knockouts may also have the potential for increased yield under high density planting.

Rice and Arabidopsis homologs of yeast CHROMOSOME TRANSMISSION FIDELITY PROTEIN 4 commonly interact with Polycomb complexes but exert divergent regulatory functions.

Both genetic and epigenetic information must be transferred from mother to daughter cells during cell division. The mechanisms through which information about chromatin states and epigenetic marks like histone 3 lysine 27 trimethylation (H3K27me3) are transferred have been characterized in animals; these processes are less well understood in plants. Here, based on characterization of a dwarf rice (Oryza sativa) mutant (dwarf-related wd40 protein 1, drw1) deficient for yeast CTF4 (CHROMOSOME TRANSMISSION FIDELITY PROTEIN 4), we discovered that CTF4 orthologs in plants use common cellular machinery yet accomplish divergent functional outcomes. Specifically, drw1 exhibited no flowering-related phenotypes (as in the putatively orthologous Arabidopsis thaliana eol1 mutant), but displayed cell cycle arrest and DNA damage responses. Mechanistically, we demonstrate that DRW1 sustains normal cell cycle progression by modulating the expression of cell cycle inhibitors KIP-RELATED PROTEIN 1 (KRP1) and KRP5, and show that these effects are mediated by DRW1 binding their promoters and increasing H3K27me3 levels. Thus, although CTF4 orthologs ENHANCER OF LHP1 1 (EOL1) in Arabidopsis and DRW1 in rice are both expressed uniquely in dividing cells, commonly interact with several Polycomb complex subunits, and promote H3K27me3 deposition, we now know that their regulatory functions diverged substantially during plant evolution. Moreover, our work experimentally illustrates specific targets of CTF4/EOL1/DRW1, their protein-proteininteraction partners, and their chromatin/epigenetic effects in plants.

Genome-wide identification of microRNAs and phased siRNAs in soybean roots under long-term salt stress.

BACKGROUND: Salinity stress, as the key limiting factor for agricultural productivity, can activate a series of molecular responses and alter gene expression in plants. Endogenous regulatory small RNAs, such as microRNAs (miRNAs) and phased siRNAs (phasiRNAs), play crucial roles during stress adaptation and prevent the injury from environmental circumstances. OBJECTIVE: To identify long-term salt stress responsive miRNAs and phasiRNAs as well as their associated genes and pathways in soybean roots. METHODS: Small RNA and degradome sequencing strategies were applied to genome widely investigate miRNAs and phasiRNAs in soybean roots under control and long-term salt stress conditions. RESULTS: In this study, stringent bioinformatic analysis led to the identification of 253 conserved and 38 novel miRNA candidates. Results of expression profiling, target and endogenous target mimics predictions provided valuable clues to their functional roles. Furthermore, 156 genes were identified to be capable of generating 21 nt and 24 nt phasiRNAs, in which 37 candidates were confirmed by degradome data for miRNA-directed cleavage. Approximately 90% of these phasiRNA loci were protein coding genes. And GO enrichment analysis pointed to "signal transduction" and "ADP binding" entries and reflected the functional roles of identified phasiRNA genes. CONCLUSION: Taken together, our findings extended the knowledge of salt responsive miRNAs and phasiRNAs in soybean roots, and provided valuable information for a better understanding of the regulatory events caused by small RNAs underlying plant adaptations to long-term salt stress.

Dissection of the Genetic Architecture of Rice Tillering using a Genome-wide Association Study.

BACKGROUND: Rice tiller number (TN) is one of the most important components associated with rice grain yield. Around one hundred rice TN genes have been identified, but dissecting the genetic architecture of rice TN variations remains difficult because of its complex trait and control by both major genes and quantitative trait loci (QTLs). RESULTS: In this study, we used a subset of the rice diversity population II (S-RDP-II), genotyped with 700,000 single nucleotide polymorphisms (SNPs), to identify the loci associated with tiller number variations (LATNs) through a genome-wide association study (GWAS). The analysis revealed that 23 LATNs are significantly associated with TN variations. Among the 23 LATNs, eight are co-localized with previously cloned TN genes, and the remaining 15 LATNs are novel. DNA sequence analysis of the 15 novel LATNs led to the identification of five candidate genes using the accessions with extreme TN phenotypes. Genetic variations in two of the genes are mainly located in the promoter regions. qRT-PCR analysis showed that the expression levels of these two genes are also closely associated with TN variations. CONCLUSIONS: We identified 15 novel LATNs that contribute significantly to the genetic variation of rice TN. Of these 15, the five identified TN-associated candidate genes will enhance our understanding of rice tillering and can be used as molecular markers for improving rice yield.

Genome-wide analysis of Dongxiang wild rice (Oryza rufipogon Griff.) to investigate lost/acquired genes during rice domestication.

BACKGROUND: It is widely accepted that cultivated rice (Oryza sativa L.) was domesticated from common wild rice (Oryza rufipogon Griff.). Compared to other studies which concentrate on rice origin, this study is to genetically elucidate the substantially phenotypic and physiological changes from wild rice to cultivated rice at the whole genome level. RESULTS: Instead of comparing two assembled genomes, this study directly compared the Dongxiang wild rice (DXWR) Illumina sequencing reads with the Nipponbare (O. sativa) complete genome without assembly of the DXWR genome. Based on the results from the comparative genomics analysis, structural variations (SVs) between DXWR and Nipponbare were determined to locate deleted genes which could have been acquired by Nipponbare during rice domestication. To overcome the limit of the SV detection, the DXWR transcriptome was also sequenced and compared with the Nipponbare transcriptome to discover the genes which could have been lost in DXWR during domestication. Both 1591 Nipponbare-acquired genes and 206 DXWR-lost transcripts were further analyzed using annotations from multiple sources. The NGS data are available in the NCBI SRA database with ID SRP070627. CONCLUSIONS: These results help better understanding the domestication from wild rice to cultivated rice at the whole genome level and provide a genomic data resource for rice genetic research or breeding. One finding confirmed transposable elements contribute greatly to the genome evolution from wild rice to cultivated rice. Another finding suggested the photophosphorylation and oxidative phosphorylation system in cultivated rice could have adapted to environmental changes simultaneously during domestication.

Dissection of the genetic architecture of rice resistance to the blast fungus Magnaporthe oryzae.

Resistance in rice cultivars to the rice blast fungus Magnaporthe oryzae is complex and is controlled by both major genes and quantitative trait loci (QTLs). We undertook a genome-wide association study (GWAS) using the rice diversity panel 1 (RDP1) that was genotyped using a high-density (700 000 single nucleotide polymorphisms) array and inoculated with five diverse M. oryzae isolates. We identified 97 loci associated with blast resistance (LABRs). Among them, 82 were new regions and 15 co-localized with known blast resistance loci. The top 72 LABRs explained up to 98% of the phenotypic variation. The candidate genes in the LABRs encode nucleotide-binding site leucine-rich repeat (NBS-LRR) resistance proteins, receptor-like protein kinases, transcription factors and defence-related proteins. Among them, LABR_64 was strongly associated with resistance to all five isolates. We analysed the function of candidate genes underlying LABR_64 using RNA interference (RNAi) technology and identified two new resistance alleles at the Pi5 locus. We demonstrate an efficient strategy for rapid allele discovery using the power of GWAS, coupled with RNAi technology, for the dissection of complex blast resistance in rice.

Development of Genome-Wide Insertion and Deletion Polymorphism Markers from Next-Generation Sequencing Data in Rice.

BACKGROUND: Next-generation sequencing technologies enable the re-sequencing of a large number of genomes and provide an unprecedented opportunity to discover numerous DNA polymorphisms throughout the genome of a species. As the second most abundant form of genetic variation, InDels, with characteristics of co-dominance, multiple alleles and high stability and density and that are easy to genotype, have received an increasing amount attention. RESULTS: In this work, a total of 2,329,544 InDels were identified in 1767 rice genomes; these InDels were dispersed across all 12 rice chromosomes, with one InDel marker found, on average, every 160.22 bp. There were 162,380 highly polymorphic InDels with a polymorphism information content (PIC) ≥ 0.5, contributing 1.81 % to the unique primer set. Of these highly polymorphic InDels, we also selected InDels with major allele differences (the size difference between the most and second most frequent alleles) ≥ 3 bp or 8 bp for primer design, which provided a more flexible choice for researchers. Finally, we experimentally validated 100 highly polymorphic InDels for accuracy and polymorphism. The PCR results showed that the accuracy of the InDel markers was 95.70 %, while the average PIC value was 0.56, with a range of 0.19 to 0.78; the average allele number was 3.02, with a range of 2 to 5. CONCLUSIONS: Our genome-wide and easily used InDel markers with high polymorphism and density in both cultivated and wild rice will undoubtedly have practical implications in rice marker-assisted breeding and will also meet the need of fine-scale genetic mapping in map-based rice gene cloning.

[Generation and identification of rice T-DNA insertional mutant lines].

By using the matured embryos of Japonica rice variety Zhonghua No. 11 as explants, rice transformation was performed by Agrobacterium-mediated co-cultivation method, resulting in 1489 independent transgenic rice plants that carry a T-DNA insertion. Genomic DNA gel-blot and PCR analyses showed that 69.8% of the total lines contain the inserted T-DNA. The flanking sequence of T-DNA in transgenic rice plants was analyzed using Tail-PCR. In addition, we have evaluated 1066 T1 transgenic lines on heading days, plant height and panicles per hill, and found different types of mutants from a number of lines.